Properties of cosmologies with dynamical pseudo Nambu— Goldstone bosons 
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We study observational constraints on cosmological models with a quintessence field in the form of 
a dynamical pseudo Nambu-Goldstone boson. After reviewing the properties of the solutions, from 
a dynamical systems phase space analysis, we consider the constraints on parameter values imposed 
by luminosity distances from the 60 type la supernovae published by Perlmutter et al., and also 
from gravitational lensing statistics of distant quasars. In the case of the type la supernovae we 
explicitly allow for the possibility of evolution of the peak luminosities of the supernovae sources, 
using simple empirical models which have been recently discussed in the literature. We find weak 
evidence to suggest that the models with supernovae evolution fit the data better in the context 
of the quintessence models in question. If source evolution is a reality then the greatest challenge 
facing these models is the tension between current value of the expansion age, Hoto, and the fraction 
of the critical energy density, Q^o, corresponding to the scalar field. Nonetheless there are ranges of 
the free parameters which fit all available cosmological data. 
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I. INTRODUCTION 

Scalar fields have played a central role in models of the 
very early universe for the past 20 years. In the past few 
years attention has turned to models in which a scalar 
field plays a dynamical role at late times, rather than sim- 
ply being frozen in as a static relic vacuum energy. Such 
models, which have been dubbed "quintessence" models 
101 , could in principle provide a dynamical solution to the 
cosmological constant problem - namely the question of 
why the magnitude of the vacuum energy at the present 
epoch is so much smaller than one might naively expect 
from particle physics models such as various supergravity 
theories. A dynamical 'solution' of the cosmological con- 
stant problem would amount to a demonstration that a 
particular dynamical evolution of the scalar quintessence 
field is a natural consequence of the cosmological field 
equations without fine-tuning of parameters, given some 
reasonable physical assumptions about the initial condi- 
tions. 

The most notable recent observational evidence which 
has driven the theoretical interest is the measurement of 
the apparent magnitude-redshift relationship using type 
la supernovae (SNe la) 1^. These results have been in- 
terpreted, in the context of a cosmological model con- 
taining pressureless dust and a cosmological constant. A, 
as evidence that the universe is undergoing accelerated 
expansion at the present epoch (see [^Q and references 
therein). The validity of this conclusion is currently open 
to some doubt, however. In particular, a recent analysis 
by Riess et al. indicates that the sample of type la 
supernovae shows a possible evolution in rise times from 
moderate (z^O.S) to large (z^l) redshifts. Although 
the statistical significance of this result has been dimin- 
ished - from the 5.8cr level § to the 1.5cr level - upon 
a more rigorous treatment of the uncertainties in the data 



it remains true that while a systematic evolution in 
the rise times of the supernovae is not conclusively ruled 
in, neither is it conclusively ruled out. 

Given that an evolution in the shape of the light curves 
of the supernovae measured in their rest frame remains a 
real possibility, it would not be surprising if the peak lu- 
minosity - which is the effective standard candle used 
- were also to evolve. Riess et al. ^ conclude that 
the type la supernovae data could conceivably be ex- 
plained entirely within the context of an open Friedmann- 
Robertson- Walker universe together with a reasonable 
astrophysical evolution model, e.g., a consequence of a 
time variation of the abundances of relevant heavy ele- 
ments in the environment of the white dwarf supernovae 
progenitors. Detailed astrophysical modelling - see, e.g., 
- should hopefully eventually resolve the issue, al- 
though at this stage the difference between our theoret- 
ical understanding and the observations remains quite 
substantial [||. 

In many recent papers it has been commonly assumed 
that the dynamical scalar field, (/>, should obey an effec- 
tive equation of state ~ wp,f, with — 1 < w < 0, at the 
present epoch, in order to obtain a cosmological acceler- 
ation, i.e., a negative deceleration parameter qo. Indeed, 
the condition that — 1 < w < is often taken as a defining 
characteristic of "quintessence" Q . The broad picture in 
this cosmological scenario is that the universe is currently 
in the early stage of an epoch of infiationary expansion. 
The motivation for this is that one could then hope to 
have a model cosmology in which observations such as the 
type la supernovae apparent magnitude-redshift relation 
could be explained by a cosmological acceleration in a 
similar fashion to models with a cosmological constant, 
but with the possibility of explaining why the magnitude 
of the vacuum energy density and the energy density in 
ordinary pressureless matter, pm, are comparable at the 
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present epoch - the so-cahed "cosmic coincidence prob- 
lem" §. 

One attractive feature of homogeneous isotropic cos- 
mological models with dynamical scalar fields is that 
many of them possess "cosmological scaling solutions" 
namely solutions which at late times have energy 
density components which depend on the cosmic scale 
factor according to p cx a~™^ and oc a"™^ simulta- 
neously, and which act as attractors in the phase space. 
If m2 < mi, which is the case, for example, for simple 
power-law potentials with inverse powers p^ , p^ ^(0) 
(/)"", or for certain power-law potentials with positive 
powers then the scalar field dominates at late times, 
producing a quintessence-dominated cosmology with ac- 
celerated expansion at late times. If mi = m2, which 
is the case for exponential potentials ||l"^-|20[|, then the 
scaling solutions are "self-tuning" jl^ - i.e., the dynam- 
ics of the scalar field follows that of the other dominant 
energy component, with a dependence oc a^"* in the 
radiation-dominated era and a dependence p^ cx a^'^ in 
the matter-dominated era. 

Even if the ultimate late-time behaviour properties of 
the solutions is not precisely "self-tuning" in the above 
sense, models such as those with inverse power law po- 
tentials can still effectively act as "tracking solutions", 
since for a wide range of initial conditions, the solutions 
rapidly converge to a common, cosmic evolutionary track 
||12| . Thus there are a number of ways in which one might 
hope to solve the "cosmic coincidence problem" , though 
in practice a degree of tuning of the parameters has been 
necessary in all models studied to date. 

In this paper, we consider a form of quintessence, an 
ultra-light pseudo Nambu-Goldstone-boson (PNGB) Q 
which is still relaxing to its vacuum state. From the view- 
point of quantum field theory PNGB models are the sim- 
plest way to have naturally ultra- low mass, spin-0 par- 
ticles and hence perhaps the most natural candidate for 
a presently-existing minimally-coupled scalar field. The 
effective potential of a PNGB field can be taken to be 
of the form 



particular, a number of constraints have been placed on 
the parameters M and / by various sets of observational 
data [0Il8|j2|J2l. Most recently, Frieman and Waga 



y(0) = M4[cos(0//) + l] 



(1) 



where the constant term is to ensure that the vacuum 
energy vanishes at the minimum of the potential. This 
potential is characterized by two mass scales, a purely 
spontaneous symmetry breaking scale / and an explicit 
symmetry breaking scale M . 

The effective PNGB mass is m^^M^/f. To obtain 
solutions with il^^l, the energy scales are essentially 
fixed to values M^IO^^ eV , interestingly close to 
the neutrino mass scale for the MSW solution to the solar 
neutrino problem, and / ~ mp^ ~ 10^^ GeV , the Planck 
scale. Since these two energy scales have values which are 
reasonable from the viewpoint of particle physics, one 
might hope to explain the coincidence that the vacuum 
energy is dynamically important at the present epoch. 

The cosmology of PNGB models has already been ex- 
tensively studied in the literature |0jl8l|2n,|^-pl . In 



have set bounds based on the SNe la data of Riess et al. 
Q (hereafter R98) on the one hand, and gravitational 
lensing surveys on the other. Comparing these bounds is 
of interest, since the Snia data have been interpreted as 
favouring a cosmological constant, whereas gravitational 
lensing data has been used to place upper bounds on A 
p5|-p7|. They therefore provide complementary tests of 
the parameter spaces of models with a non-trivial vac- 
uum energy. 

In this paper it is our intention to critically study these 
bounds. Firstly, we will consider how the bounds are af- 
fected by the initial value of the scalar field at the begin- 
ning of the matter-dominated epoch. Secondly, we wish 
to investigate how such bounds might be affected in the 
case of the PNGB model if the observed apparent faint- 
ness of type la supernovae is at least partly due to an 
intrinsic evolution of the sources over cosmological time 
scales, which in view of the results of Q would appear 
to be a very real possibility. The reason for focusing on 
the PNGB models in such an investigation is suggested 
by the fact that whereas many quintessence models have 
been singled out in the literature, perhaps somewhat ar- 
tificially, simply because they have the property of yield- 
ing an accelerated expansion, many different possibilities 
arise in the PNGB case. Indeed, at very late times, the 
apparent magnitude-redshift relation for PNGB models 
ultimately coincides with that of the Einstein-de Sitter 
model, even though the density of ordinary matter can be 
low in the PNGB cosmologies. The requirement that the 
faintness of type la supernovae is entirely due to their 
cosmological distances places rather strong restrictions 
on the values of the parameters M and / |2^] , because it 
requires us to exist at an epoch of the PNGB cosmologies 
which is still quite far removed from our ultimate destiny. 
If these restrictions are relaxed because of evolutionary 
effects, then it is quite plausible that other regions of the 
parameter space of the PNGB models become viable al- 
ternatives. Since PNGB cosmologies could therefore still 
solve the "missing energy problem" , even if the evidence 
for a cosmological acceleration proves to be ephemeral, 
we believe it is important to investigate this possibility 
quantitatively. 

We will begin the paper with a qualitative analysis of 
the solutions, to provide some general insights which will 
help to guide our quantitative discussion. Although these 
properties are no doubt already known, to the best of our 
knowledge an analysis of the phase space of the solutions 
has never been presented in the literature. Having com- 
pleted this analysis in Sec. II we will go on to discuss 
a number of issues relating to numerical integration in 
Sec. Ill, and relate the properties of the solutions found 
numerically to the exact analysis of Sec. II. In Sec. IV we 
present the main analysis of the constraints imposed on 
the (M, /) parameter space, allowing for the possibility 
of evolution of peak luminosities in the type la supernova 
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sources. Bounds from gravitational lensing statistics are 
updated in Sec. V, and the implications of our results are 
discussed at greater length in Sec. VI. 



II. PHASE-SPACE ANALYSIS 

We will begin by performing an analysis of the differ- 
ential equations governing the cosmological evolution in 
a manner similar to p revious studies in inflationary and 
quintessential models 
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The classical action for gravity coupled to a scalar field 
has the form 



S = 



c 



(2) 



where k is the Planck constant, R is the Ricci scalar, 
g = detg^u, and L is the Lagrangian density of non- 
relativistic matter and radiation. For simplicity, we as- 
sume (j) is minimally coupled to the curvature, and we 
work in units in which ?i = c = 1 . 

Consider a spatially-flat Friedmann-Robertson- Walker 
(FRW) universe containing a fluid with barotropic equa- 
tion of state P~f = (7 — i)p~f, where 7 is a constant, 
< 7 < 2, such as radiation (7 = 4/3) or dust (7 = 1). 
There is a self-interacting scalar field with the PNGB po- 
tential energy density (m) evolving in this universe. The 
total energy density of this homogeneous scalar field is 
P4> — 0^/2 + V{cl)). The governing equations are given by 



H = 
P7 = 



p. 



-3H{p^ + P^) , 



subject to the Friedmann constraint 



H' = 



1 



Pi 



V 



(3) 
(4) 
(5) 

(6) 



where = SttG, H — a/ a is the Hubble parameter, and 
an overdot denotes ordinary differentiation with respect 
to time t. 

We may rewrite the Friedmann constraint as 



where 



K^Pj 

3IP ' 



1 



3iJ2 V 2 



(7) 

(8) 
(9) 



are the ratios of the energy densities of the barotropic 
matter and the quintessence field as fraction of the criti- 
cal density respectively. 



In contrast to the case of the cosmologies with an ex- 
ponential potential [|l4|Jl5| , [l9[| where the dynamics can be 
reduced to a 2-dimensional autonomous phase plane, for 
the system - (||) the simplest phase space appears to 
be 3-dimensional in the full four-dimensional phase space. 

There are two alternative choices of variables which are 
useful to describe the dynamics, which we will discuss in 
turn. 



A. Field Variables 

The first choice is to simply use the Hubble parameter, 
H, the scalar field and its first derivative as the elemen- 
tary variables. These are of course simply the variable u, 
w, w of Frieman and Waga up to an overall scaling. 
By defining 



/ = 



/ 



Kcf) 



J = K,(j) ; 



we therefore obtain the system 



2 



7m 



(cos / 



J 



-iHJ 



sin/ , 



(10) 

(11) 
(12) 

(13) 



where for notational simplicity we define rv? = k^M"*, 
and J- = nf . so that T is dimensionless, while m has di- 
mensions inverse time. The constraint equation becomes 



P-y 



1 



m^{cosI + 1) J 



(14) 



From Eq. it follows that py = if = 0. 
Therefore trajectories do not cross the 2-dimensional 
p-y — surface, which is a hyperboloid in the variables 
H, cos(//2), and J. Physical trajectories with p-^ > 
are forced to lie within the volume of the H, I, J phase 
space bounded by the p-y = surface. 

The only critical points of the system ([Tl])-(p^ at finite 
values of H, I, J occur at 

1. Ci± a.t H ^ ±Hi, / = mod 27r, J = 0; and 

2. C2 at iJ = 0, / = TT mod 27r, J = 0, 
where 



Hi 



(15) 



Both of these points in fact lie on the = surface. Fur- 
thermore, this surface intersects the H — plane only at 
the isolated points C2. The H > and H < subspaces 
are thus physically distinct, and the H < subspace 
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simply correspond to the time-reversal of the H > sub- 
space. Therefore we can take H > without loss of 
generality. 

The pattern of trajectories close to the p-y = surface 
can be ascertained by continuity to the p-y ~ solutions, 
even though the latter are not physical. The — 
subspace is obtained, for example, by regarding (jlj) as 
a quadratic equation for H, and using the solution to 
eliminate H, thereby obtainin g a 2-dimensional system 
for / and J given by (^2|) and (p^). 

We plot the resulting H > pattern of trajectories in 
Fig. |l| for values of / e [0, 2tt). Since the potential, V{(f>), 
is periodic the same pattern of trajectories repeats itself 
as we extend / to ±oo, with trajectories crossing from 
one "cell" to another at the cell boundaries. 




FIG. 1. The projection of the trajectories within the 
= subspace on the I - J plane for values of 7 £ [0, 2tt). 
Within this subspace, Ci+ is a saddle point and C2 is a stable 
spiral. 

The trajectories which occupy the lower half of Fig. |l| 
are obtained from those in the upper half by the sym- 
metry / — > 27r — J, J — > — J of the differential equations 
( |lT| ) - (|lj). Physically this simply corresponds to the 
scalar field rolling from the maximum to the right of the 
minimum as opposed to the one to the left. 

An analysis of small perturbations about the critical 
points Ci and C2 yield eigenvalues 



1. A = -\/6TO7,-f {y6± v'eTIT^ 



at Ci 



2. A = 0, 



I i7n 



at C2 



Thus Ci+ attracts a two-dimensional bunch of trajec- 
tories but is a saddle point with respect to trajectories 
lying in the p-y = surface, as is evident from Fig. ||. 
The 2-dimensional bunch of trajectories which approach 
Ci-(- are found to correspond to inflationary solution with 
a (X exp(-y/2/3mi) as t ^ 00 and (j) ~^ const = 2mTJ-, 



n G Z. The possible role of scalar fields with PNGB po- 
tentials in driving an inflationary expansion of the early 
universe has been discussed in [ psj . 

The point C2 is a degenerate case, in particular with 
regard to perturbations orthogonal to the p-y — surface, 
(i.e. into the surface p^ > region), for which the eigen- 
value is zero. It is a centre with respect to the trajecto- 
ries lying in the p-y = surface, and when perturbations 
of higher order are considered it becomes a stable spiral 
point in the p^ = surface as can be seen in Fig. |^. Since 
there is a degeneracy, however, an alternative choice of 
phase space variables is desirable. We will defer a discus- 
sion of the late time behaviour of the solution near C2 to 



Sec. II B 



The points Ci± correspond to models with a scalar 
field sitting at the maximum of the potential, whereas 
C2 corresponds to the scalar field sitting at the bottom 
of the potential well. The separatrices in Fig. ^ which 
join Ci± to C2 correspond to the field rolling from the 
maximum to the minimum. It would appear from Fig. |l| 
trajectories which spiral into C2 become arbitrarily close 
to the separatrix at late times. 

The separatrices which join the points Ci± to points 
at infinity correspond to solutions for which the scalar 
field reaches the top of the potential hill as t — > ±00, 
(e.g., the rightmost trajectory in Fig. |l|). Finally, there 
are also straight line separatrices parallel to the 77- axis 
at each of the points Ci±, extending from 77 = ±77i 
to infinity, which represent solutions with a static scalar 
field sitting on top of the potential hill. 

To examine the critical points at infinity it is conve- 
nient to transform to spherical polar coordinates r, 9, 
and (p by defining 



H — r cos 9 , 
I = r sin 9 sin (p , 
J = r sin 9 cos 



(16) 
(17) 
(18) 



and to bring the sphere at infinity to a finite distance 
from the origin by the transformation r = p/{l — p), 

o< p< 1 11. 

Although the trajectories on the sphere at infinity do 
not represent physical cosmologies, it is useful to plot 
them since the form of the trajectories which lie just 
within the sphere will be similar. On the sphere p = 1 
we find 



d9_ 



sm f cos 



3 sin 



7 



(tan^ 9 cos^ </> - 6) 



— — - cos 9 sm 20, 
d( 2 ^ 



(19) 
(20) 



where ^ is a new time coordinate defined by d^ = rdt. 
The resulting integral curves are plotted in Fig. |. By 
( p^ ) the projection of the physical region py > onto 
the sphere at infinity leads to the condition 
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cot^ 6 > - cos^ I 
6 



(21) 



Values of 6 and (j) which violate this inequality lie in the 
shaded region. 

The critical points on the sphere at infinity are 

1. A^j^, A^: four points at 

{e, 0) e {(± tan-i Ve, 0), (± tan-i ^6, tt)} or 
H = ±oo, I/H = 0, and J/H = ±V6. 

2. B^': two points at 

H = ±oo, I/H = 0, and J/H = 0. 

Since the projection onto spherical polar coordi- 
nates (jl^) - ( p^ is degenerate at the north and 
south poles 6 — 0,TT, these points are excluded from 
the chart {9, (j)) but can be included using an alter- 
native hemispherical projection (see Fig. ^). 

3. CJ°2- two points at 
(0,0)G{(7r/2,V2),(7r/2,3V2)}or 
/ = ±oo, H/I = 0, and J/I = 0. 

TABLE I. The critical points on the sphere at infinity and 
their eigenvalues. 



Critical points I 



Eigenvalues (with degeneracies) 



Coc 
1,2 



±3 (2), ±3(2-7) 



±^ (2), t|(2-7) 



0(3) 



An analysis of small perturbations shows that the 
points A^ 2+ ^-re repellors in all directions of the phase 
space (cf., Table |), while A^j^ 2- 

are attractors. This 
therefore represents the most 'typical' early behaviour 
of solutions. For H > 0, 2+ correspond to the 
limit t ^ 0. We find that or a cx t^^^, while 

±y'2/31nt, for these solutions. The points A^j^ 2- 
with H < represent the time-reversed solutions. 

The point B!^ (B'^) repels (attracts) a 2-dimensional 
bunch of trajectories travelling to (from) finite values of 
H, I, J, but is a saddle point with respect to directions on 
the sphere at infinity. The points are found to correspond 



tot^O with iJ- 3^ or a cx t2/37 ^hile cx t", n > 0. 

The points the projection of the points Ci± 

and C2 into the sphere at infinity. The degenerate eigen- 
values simply refiect the degeneracy of the projection. 

Bjp acts as a repellor for trajectories with ~ 0, (/) ~ 
const ast ^ 0. As shown in Fig. |[ trajectories are driven 
towards B"^ before they reach CJ°. This is consistent 
with the property that when H is large {3H > m^), the 
field evolution is over-damped by the expansion, and the 
field is eflFectively frozen to its initial value ((/) 0). 



B. Energy Density Variables 

In view of the eigenvalue degeneracy encountered 
above, we can alternatively choose to represent the sys- 
tem by the variables H , x, and y, where 



y 



K(j) J 

V6H ~ VGH ' 
kVV \/2tocos(//2) 



ViH 



VSH 



(22) 
(23) 



These are the same variables used by Copeland, Liddle 
and Wands [|l9| for the model with an exponential po- 
tential. As above, we will consider H > Q only. The field 
equations then take the form 



H = — , 



(24) 



my j 3y^H'^ 3 



2m^ '2^-(A^-2), (25) 



mx I 3y'^H^ 3 
y = T — \/ 1 + -^Hyfj. , 



where 



fi{x,y)=-f{l-y^) + i2-j)x^ 
We note that in these variables 



y 



(26) 



(27) 



(28) 



which is why we have adopted the terminology "energy 
density variables" . The physical region of the phase space 
will be constrained to lie within the cylinder x'^ + y^ < 1 
since il^ < 1- In the case of the exponential potential 
analysed in Ref. one of the differential equations 

decoupled, and the dynamics was eflFectively described 
by a phase plane with trajectories bounded by the cir- 
cle x"^ + y"^ = 1. In the present case, however, no such 
simplification arises. 

The physical region of phase space is further restricted 
by the requirement that < ^m? /H^ = Hi^ / H^, which 
is equivalent to cos^ (^/2) < 1 in terms of the field vari- 
ables. For values H > Hi each H = const slice of the 
cylinder a:^ -I- < 1 is cut off in the y-direction above 
and below the y = ±Hi/H lines. Thus the "fundamental 
cell" of the phase space can be considered to be a cylinder 
for Q < H < Hi, capped by a horn for H > Hi, which 
tapers off to a line segment — 1 < a; < 1 on the x-axis as 

— > 00 (see Fig. In fact, the phase space consists 
of an infinite number of copies of the fundamental cell of 
Fig. 1^ as a result of the periodic structure of the potential. 
These cells, C„, can be labelled by an integer, n, with the 
variable / lying in the range 2mT < / < 2(n -I- l)7r for 
each n. For cells with even n the dynamics is described 
by (p^-(27) with the upper sign in (^) and (26), while 
for odd n one must take the lower sign. 
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FIG. 2. The projection of trajectories within the sphere at infinity on the 4> - 9 plane. The unphysical region is shaded. 




FIG. 3. The projection of trajectories within the / > and / < hemispheres at infinity on the H — J plane. The 
unphysical region is shaded. 
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FIG. 4. The "fundamental cell" of the phase space in 
terms of the energy density variables. On the = Hj/H^ 
planes (The shaded planes), only one trajectory is possible as 
shown. It corresponds to the scalar field lying on top of the 
potential hill all along. On the H = plane, the trajectories 
are concentric circles with centre at the origin. 



y — A{t) cos 



mt 



(30) 



(31) 



where the amplitude is governed by the equation 
The nature of C2 is now found to depend on 7: 



1. If 7 < 1 we find that over a cycle the average value 
of the r.h.s. of ( |3l] ) is negative and A? decreases 
so that C2 is an asymptotically stable spiral. Fur- 
thermore to leading order H ~ as t 00 or 

a cx t^^"^^ and A{i) is given by 



A(t) = Bt(^-i)/^exp 



1 

7 Jt 



dt' f mt' 
— cos — — 
t' V ^ 



(32) 



with B constant. The late-time attractor has 



and fl^ 
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2. If 7 > 1 then over a cycle the average value of the 
r.h.s. of ( |3l] ) is positive and A^ increases until it 
reaches a limit cycle A^ — 1, i.e. + ^ 1 or 



1, r^T, = 0. In this case 



Within the horn portion of each phase space cell 
the motion of most trajectories is roughly circular in 
H = const slices, in a clockwise sense in even cells, 
and anti-clockwise in odd cells C2n+i- However, tra- 
jectories can cross from one cell to another along the 
y = ±Hi/H boundaries of their horns, which correspond 
to the surfaces / = 2mT in terms of the field variables. 
For even cells, C2n, trajectories join the cell C2n along 
the y = Hi/H surface and the cell C2n+i along the 
y = —Hi/H surface. Below the H = Hi plane solu- 
tions cannot cross from one cell to another, but remain 
confined within the cylinder + y"^ < 1. 

When H ^0 we see from (|2|) - (|^ that x = my/T 
and y = —mxjT^ so that trajectories which lie in the 
H = surface are purely concentric circles. Since d = 
in the H — plane, these do not represent physically 
interesting cosmologies, but by continuity the behaviour 
of the trajectories just above the plane will be of a spiral 
nature. 

The origin H = x = y = is in fact the critical point 
corresponding to C2. The nature of the critical point is 
altered by the change of variables, however. In particu- 
lar, whereas the eigenvalues for linear perturbations are 
unchanged, when higher order corrections are considered 
the point is no longer always an asymptotically stable 
spiral as was the case in Fig. 0. 

Asymptotically as i — > 00 we have 



A{t) = 1 - Bi2(i-7) 



exp 



X — A{t) sin 



mt 



(29) 



r „ r dt' 






^mt'J 







(33) 



3. If 7 — 1, which corresponds physically to an or- 
dinary matter-dominated universe, then an inter- 
mediate situation obtains. Essentially any of the 
concentric circles in the H — plane of Fig. || can 
be approached asymptotically giving a universe for 
which — > ai and I — ai where ai is a 

constant in the range < ai < 1, which depends 
on the initial conditions and the parameters m and 



We observe that for all values of 7, 



1 



(34) 



at late times. Since cx a^^'^ , the three different late 
time behaviours can thus be understood as a consequence 
of the scalar field either decreasing more rapidly than 
the barotropic fluid (7 < 1), less rapidly (7 > 1), or at 
the same rate (7 = 1). The scalar field thus eventually 
dominates if 7 > 1, while the barotropic fluid dominates 
if 7 < 1. In the interesting critical case of dust filled 
models (7 = 1) both the scalar field and ordinary matter 
are of cosmologically significant density at late times. 

One simplification that is often made in studying 
quintessence models is to assume that at late times the 
quintessence field obeys an equation of state 



7 



^ = (7^-1)/ 



(35) 



with 70 effectively constant. While such an assumption is 
justified in the case of models with a slowly varying scalar 
field, it docs not apply in the present case. In particular, 
since the effective barotropic index of the scalar field is 



70 



2x^ 



(36) 



we see that 7^ ~ sin^ {mt/T) at late times, so that it 
remains truly variable, varying from to 2 over each 
cycle. It follows from ([25|)-([2^) that the scalar energy 
density parameter obeys the equation 



Q.^ = iHVt^ (1 - n^) (7 - 70) 



(37) 



so that — > as t — > cx), which accords with the late 
time properties of the solutions observed above. 



III. NUMERICAL INTEGRATION 

We will now consider the extent to which models based 
on the PNGB potential are constrained by the latest ob- 
servational evidence. This work extends the studies pre- 
viously undertaken by various authors p7| , p^ , p3|j24| . In 
the most recent analysis, Frieman and Waga have com- 
pared the constraints imposed by the high redshift su- 
pernovae luminosity distance on the one hand and grav- 
itational lensing bounds on the other. The two measures 
provide tests which are potentially in opposition. Here 
we will perform a similar analysis for the PNGB models, 
but also taking into account the possibility of luminos- 
ity evolution which has not been considered in previous 
studies 24 1 . We are therefore considering the model 
of the previous section with 7 = 1. 

To proceed it is necessary to integrate the equations 
numerically. To do this we introduce the dimensionless 
variables 



Ho Ho 

K(t) 

W = 1 = — — . 

J" 



(38) 

(39) 
(40) 



where fimo is the fractional energy density of matter at 
the present epoch, to. The dynamical system then be- 
comes 



.H m2 . 



V = —V 



Ho ' :fhi 

H 

H~o ' 



w — — , 



(41) 

(42) 
(43) 



where the Hubble parameter is defined implicitly accord- 
ing to 



H 

H^ 



1 2 



3i72 



{cosw + 1) 



(44) 



and prime denotes a derivative with respect to the di- 
mensionless time parameter r = Hot. Only the variable 
V differs from those used by Frieman and Waga ||l^ . Our 
reason for making the choice (39) is that it allows us to 
integrate the Friedman equation directly rather than a 
second order equation ( [lT| ) which follows from the other 
equations by virtue of the Bianchi identity. This may 
possibly lead to better numerical stability since it is not 
necessary to implement the Friedmann constraint sepa- 
rately. 

We begin the integration at initial values of u, v, and 
w chosen to correspond to initial conditions expected in 
the early matter-dominated era. The integration pro- 
ceeds then until the r.h.s. of (^J) is equal to 1, thereby 
determining the value of the present epoch, to, to be the 
time at which H = Ho- We are then also able to deter- 
mine flmo, since according to (39) 



mO 



vHto) 



(45) 



The choice of appropriate initial conditions has been 
previously discussed p^ , p^ , p^ . In particular, since the 
Hubble parameter is large at early times, it effectively 
acts as a damping term in Eq. (|^), driving the scalar 
field to a state with (p = initially, i.e. u — 0. We 
take V = 1101 initially, which in view of relation ( |39| ) 
and the fact that 0.1 — 1 corresponds to the early 

matter dominated era 1100 < z < 3000. Results of the 
integration do not change significantly if v is altered to 
values within the same order of magnitude. 

The initial value of the scalar field variable, Wi = 
w{ti), can lead to some variability in predicted cosmo- 
logical parameters at the present epoch. A few differ- 
ent values of w(ti) have been considered by different au- 
thors | p^ , p7| , p^ . However, the only systematic studies of 
bounds in the Af ,/ parameter space have been performed 
| p^ , p4[ for one particular initial value, Wi = 1.5. One 
must bear in mind that such bounds are also dependent 
on Wi, the value of which is not greatly restricted. Given 
that we are starting with u(ti) ~ 0, so that the kinetic 
energy of the scalar field is initially negligible, the only 
physical restriction on the value of Wi comes from the re- 
quirement that the scalar field should be sufficiently far 
from the minimum of the potential, V{(j)), that ri0(ii) is 
small. Thus will ensure that Wi is consistent with a scalar 
field that has emerged from the radiation dominated era 
with fl^ sufficiently small that is consistent with bounds 
set by primordial nucleosynthesis, and by structure for- 
mation models. This still leaves considerable latitude for 
the choice of Wi, however. 

In Figs. H and ^ we display contour plots of il^o and 
Hoto in the M,f parameter space for two values Wi — 
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Wi = 1.5 
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(b) M (lO"V/^eV) 



FIG. 5. Contours of ^4,0 in the M,f parameter space for two choices of initial values: (a) Wi = 1.5; (b) Wi = 0.5 
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(a) M (lO^V/^'eV) 




FIG. 7. Contours of go in the MJ parameter space for two choices of initial values: (a) Wi = 1.5; (b) Wi = 0.5. Values 
qo < 0, corresponding to a universe whose expansion is accelerating at the present epoch, are shaded. 
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1.5 and Wi = 0.5. Similar figures have been given by 
Frieman and Waga |l^ in the Wi = 1.5 case, although 
our resolution is somewhat better. As Wi decreases the 
contour plots do not change significantly in terms of their 
overall features, but contours with equivalent values shift 
to lower values of the / parameter. For example, for 
large values of M the ft^o — 0.7 contour lies at a value 
/ ~ 2.05 X 10^* GeV if Wi = 1.5, while the same contour 
lies at / ~ 0.94 x lO^* GeV if w, = 0.5. 

The other principal feature of the plots |^ and ^, which 
was not commented on in Ref . p^ ] , is the wave-like prop- 
erties of the contours at larger values of M. These fea- 
tures can be readily understood by considering the corre- 
sponding plots of the deceleration parameter, qq, which 
is defined in terms oiu{to), v(to), and t«(io) by 



3 1 
= i (1 + Sn^o - 6z/o') . 



(46) 



We display contour plots of go in the M,f parameter 
space in Fig. |^. Essentially, as M increases for roughly 
fixed / for sufficiently large M, the value of qq oscil- 
lates over negative and positive values, from a minimum 
of go = 5 (1 ~ Sfl^o) to a maximum go = ^ (1 + 3il</,o) 
about a mean of go — 0.5. This corresponds to the scalar 
field, (j), having undergone more and more oscillations by 
the time of the present epoch. The minimum value of go 
is attained when (j) — instantaneously, while the max- 
imum value of go is attained when is instantaneously 
passing through the minimum of its potential. 

For smaller values of M to the left of the plots the 
scalar field has only relatively recently become dynam- 
ical, whereas for larger values of M, the scalar can al- 
ready have undergone several oscillations by the time of 
the present epoch, particularly if / is small. This varia- 
tion can be understood in terms of the asymptotic period 
of oscillation of solutions which approach Ci, which by 
Eqs. @, {M is 



(47) 



The period ta is shorter for larger M, or for smaller /. 
Since the final / values plotted in the Wj = 0.5 case 
are a factor of 2 smaller than the Wi = 1.5 case, this 
also explains why points with the same value of M have 
undergone more oscillations up to the present epoch for 
the smaller value of Wi . 

The value of i?o^o oscillates as M increases for roughly 
fixed /, according to whether the universe has been ac- 
celerating or decelerating in the most recent past, with 
more rapid variation for parameter values with shorter 
asymptotic periods, ta- 

The wiggles in the V,^o contours are a residual effect 
of the oscillation of the scalar field as 17^ settles down 
to a constant value according to (^T\). The variation in 
the value of Q^q with the value of Wi can be understood 
from the fact that for smaller values of w, the scalar field 



begins its evolution in the matter-dominated era closer to 
the critical points, Ci+, corresponding to the maximum 
of the potential, V{(l)). For fixed M and / the period of 
quasi-infiationary expansion is therefore longer, and the 
present value of fl^o larger. 



IV. CONSTRAINTS FROM HIGH-REDSHIFT 
TYPE lA SUPERNOVAE 

Empirical calibration of the light curve - luminosity re- 
lationship of type la supernovae provides absolute mag- 
nitudes that can be used as distance indicators. Since 
the luminosity distance of a light source, cIl, is defined 
by 



Hodi 



(z + 1) / (z + l)dT, 



(48) 



with T ~ Hot, it is convenient to define an additional 
variable 



r = — I V dr 

'TO 



(49) 



which is proportional to the comoving coordinate dis- 
tance 



(50) 



' comoving — „ U \ /j- \ ' 

We then adjoin a differential equation 



(51) 



to the differential equations (|3^)-(40) when performing 
the numerical integration. In terms of r and v, the lumi- 
nosity distance is then determined according to 



HodL 



2/3 



(r-rito)) , 



(52) 



and can be used in the appropriate distance modulus to 
compare with the supernovae data. 

Frieman and Waga have recently considered con- 
straints on PNGB models from supernovae data (without 
source evolution) using the data set of Riess et al. [|| . We 
will perform a similar analysis, but we will make use of 
the largest available data set, namely the 60 supernovae 
published by Perlmutter et al. [|| (hereafter P98). Of 
these, 18 low redshift SNe la were discovered and mea- 
sured in the Calan-Tololo survey |^l| , and the Supernova 
Cosmology Project discovered 42 new SNe la at redshifts 
between 0.17 and 0.83. The peak magnitudes of the su- 
pernovae are corrected using the "stretch factor" light 
curve fitting method |^,^. This method is based on 
fitting a time-stretched version of a single standard tem- 
plate to the observed light curve. The stretch factor is 
then used to estimate the absolute magnitude. 
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FIG. 8. Confidence limits on M J parameter values, with uii = 1.5 and no evolution of sources, for the 60 supernovae la in 
the P98 dataset. Parameter values excluded at the 95.4% level are darkly shaded, while those excluded at the 68.3% level are 
lightly shaded. 



A. Models without evolution 

We use the stretch-luminosity corrected effective B- 
band peak magnitude in Table 1 of P98 as the absolute 
magnitude and denote it as tu^b ■ The distant modulus 
of a supernova is defined as rrig + A^O: where A^o is the 
fiducial absolute magnitude which has not been given in 
the literatures. A fitting method to obtain A^Oj using 
only the 18 Calan/Tololo supernovae, can be found in 
[^,|^. Perlmutter et al. |^ calculate the confidence re- 
gions by simply integrating over M.q. In this paper, we 
will perform an analytic marginalization over A^p and 
obtain the marginal likelihood. 

Similarly to the statistic used by Riess et al. Q , we 
define the quadratic form 



60 

E 

i=l 



Mo- Hi 



(53) 



where ^i{zi; Hq,M, /) is the predicted distance modulus 
for model parameters (M, /) , and 



5 loge 



The predicted distance modulus is 

= b\ogdL{z^) + 25 , 



(54) 



(55) 



if the luminosity distance, d^, as defined by ( [48| ) or (52), 
is given in units of Megaparsecs. 

In order to perform analytic marginalization over Hq 
as well as over AIq we separate out the Hq and A^o de- 
pendence from into a quantity, which we define by 



- Mo = gi 



(56) 



where gi{zi] M, f) depends implicitly only on M and /. 
We then follow the statistical procedures adopted by 
Drell, Loredo and Wasserman |35 3^ and marginalize 
over v using a flat prior that is bounded over some range 
Ai/. 

The marginal likelihood is 



1 



where 



Q{MJ) = 



dv e-xV2 



-Q/2 



2=1 



(57) 



(58) 



and i>(Af, /) is the best-fit value of given M and /, 
with conditional uncertainty s. We identify 1/s^ as the 
coefficient of i^^ in (p3|) , 



1 



60 ^ 



and i) as times the coefficient of —v, viz. 

9t 



(59) 



(60) 



The quadratic form (|5q ) is what one would obtain 
by calculating the "maximum likelihood" for M and 
/. Since s is independent of M and /, it follows from 
Eq. (57) that the marginal likelihood is proportional to 
the maximum likelihood in this case. 



i=l 



13 



Fig. H shows the 68.3% and 95.4% joint credible regions 
for M and /, and is a direct analogue of Fig. 1 of Ref. , 
where a similar analysis was performed on 37 supernovae 
given in R98. The position of the region of parameters 
which are included at both the 68.3% and 95.4% confi- 
dence levels is broadly similar to that obtained from the 
R98 data [||| . Although Frieman and Waga |||| did not 
include the parameter region M > 0.004/i eV , we have 
redone their analysis on the R98 data and find that the 
parameter values to the right of Fig. || which are admit- 
ted at the 95.4% level but excluded at the 68.3% level for 
the P98 dataset, (labelled region II in Fig. H), are in fact 
excluded also at the 95.4% confidence level if the 37 su- 
pernovae of the R98 dataset are used. It is possible that 
this discrepancy has its origin in the different techniques 
used by Riess et al. to determine the distance moduli. 
Possible systematic discrepancies in the "stretch factor" 
method of P98 versus the "multi-colour light curve" and 
"template fitting methods" of R98 have been discussed 
in some detail in Ref. |35[ . 

The importance of the 2a included parameter region 
to right of Fig. p diminishes, however, if one compares it 
with Figs. ^ and 6[ since it largely corresponds to param- 
eter values with fi^o ^ 0.9, which can be discounted by 
dynamical measurements of 57mO, where J7mO + f^0o = 1- 
Furthermore, the few allowed values below the fi^o = 0.9 
contour in this part of the parameter space have unac- 
ceptably small values for the age of the Universe, HqIq. 

The parameter region 0.002/i eV < M < 0.003/i eV 
which from Fig. || is admitted at both the 68.3% and 
95.4% levels, (labelled region I), by contrast, corresponds 
to acceptable values of both Jl^o and HqIq. Comparing 
with Fig. |7, we see that this region has —0.1 <qo< —0.6, 
corresponding to a Universe with a scalar field still in an 
early stage of rolling down the potential V{(j)) . The label I 
is thus indicative of the fact that the scalar field is rolling 
down the potential for the first time (from left to right), 
while in region II the scalar field is rolling down the po- 
tential for the second time (from right to left). In region 
II go is positive - however, it corresponds to parameter 
values for which there would have been a cosmological ac- 
celeration at modest redshifts in the past, e.g., at z ~ 0.2, 
well within the range of the current supernovae dataset. 

Since conclusions regarding statistically preferred re- 
gions of the parameter space can change if evolution of 
the sources occurs, we think it is important that this 
possibility is also examined, as we will now do. 



precisely how the source peak luminosities vary with red- 
shift, one approach is to assume some particular empir- 
ical form for the source evolution, and to examine the 
consequences. Such an analysis has been recently pre- 
formed by Drell, Loredo and Wasserman in the case 
of Friedmann-Lemaitre models with constant vacuum en- 
ergy. We will undertake an equivalent analysis for the 
case of PNGB quintessential cosmologies. 

Following Drell, Loredo and Wasserman we will 
assume that the intrinsic luminosities of SNe la scale as 
a power of 1-1- 2 as a result of evolution. This model intro- 
duces a continuous magnitude shift of the form /31n(l-|-z) 
to the SNe la sample. Eq. (|53|) then becomes 



X'(M,/) 



{rn\ 



5,-i/-/31n(l + z,))' 



(61) 



The parameter (3 will be assumed to have a Gaussian 
prior distribution with mean (3q and standard deviation 
h. Physically the parameter /3o represents a redshift- 
dependent evolution of the peak luminosity of the su- 
pernovae sources, which might be expected to arise as 
a result of the chemical evolution of the environment of 
the supernovae progenitors as abundances of heavier ele- 
ments increase with cosmic time. Ultimately, one should 
hope to account for this evolution with astrophysical 
modelling of the supernovae explosions 0|. The param- 
eter h would then account for a local variability in the 
supernovae environments between regions of individual 
galaxies at the same redshift which are richer or poorer 
in metals, or with progenitor populations of different ages 
and masses etc. 

We now have two parameters to marginalize over, v 
and (3. As in the case of models with no evolution, we will 
marginalize over v using a flat prior. We use a Gaussian 
prior for (3 with mean /3o and standard deviation 6, so 
that 



1 



6V27r 



where 



(62) 



(63) 



The marginal likelihood is calculated by multiplying 
the prior (^) by the likelihood resulting from Eq. (|5j 
and integrating over /?. The resulting likelihood is 



B. Models with evolution 



In view of the recent results of Refs. [^|Jg] , supernovae la 
may exhibit some evolutionary behaviour, at least as far 
as their rise times are concerned. A possible evolution in 
the peak luminosity is therefore a possibility which must 
be seriously investigated. 

In the absence of a detailed physical model to explain 



C{M, /) 



where 
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(65) 



14 



2 4 6 8 10 

(a) ^„ = 6 = M (10-V/^eV) 




(c) ^„ = (i\ b = 0.5 M (10-V/^eV) 



FIG. 9. Confidence limits on MJ parameter values in the best-fit /3o = 0.414 slice of the {M,f,l3o) parameter space, with 
Wi = 1.5, for tlic 60 supernovae la in the P98 dataset. Parameter values excluded at the 95.4% level are darkly shaded, while 
those excluded at the 68.3% level are lightly shaded. For reference, contours of Q^o and Hoto are superposed as dashed and 
dotted lines respectively. 
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(b) ^0 = 0, 6 = 0.25 M (10-V/^eV) 




FIG. 10. Confidence limits on MJ parameter values in the /3o = slice of the {M,f,Po) parameter space relative to a 
best-fit value /3* ~ 0.414, for the 60 supernovae la in the P98 dataset, with Wi = 1.5. Parameter values excluded at the 95.4% 
level are darkly shaded, while those excluded at the 68.3% level are lightly shaded. For reference, contours of $7,^0 and Hoto 
are superposed as dashed and dotted lines respectively. 
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(a) = 0.435, 6 = M (iQ-V/^eV) 




2 4 6 8 10 

(b) = 0.435, 6 = 0.25 M (lO-V/'eV) 




7 4 6 8 10 

(c) = 0.435, 5 = 0.5 M {iQ-'h'^'eV) 

FIG. 11. Confidence limits on M,f parameter values in the /3o = slice of the (M, /, /3o) parameter space relative to a 
best-fit value /3* ~ 0.435, for the 60 supernovae la in the P98 dataset, with Wi = 0.2. Parameter values excluded at the 95.4% 
level are darkly shaded, while those excluded at the 68.3% level are lightly shaded. For reference, contours of Q^o and Hoto 
are superposed as dashed and dotted lines respectively. 
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h^{z,- M, f) = - 5^ - A ln(l + z,) , (66) 



The conditional best-fit of /3' is given by 



(67) 



60 



/3(M, f)=a^Y. 



i=l 



ln(l + z,) 



(68) 
(69) 



and the P' uncertainity, a, is given by 



1 1 2r2 



E 



ln'(l + 2^) 



(70) 



Although (7 is independent of M and /, the marginal 
likelihood is no longer proportional to the profile likeli- 
hood because Q is now given by ( |65| ) rather than by the 
chi-square type statistic (|5^). 

We have performed a detailed numerical analysis on 
the P98 dataset, varying /3o, h and wi. As a result we 
find a best-fit value of /Jo = /3* — 0.414, for Wi = 1.5 in 
the PNGB models. This would correspond to supernovae 
being intrinsically dimmer by 0.17 magnitudes at a red- 
shift of z = 0.5, which is an effect of the typical order of 
magnitude being addressed in current attempts to better 
model the supernova explosions 0|. Furthermore, we find 
that inclusion of a non-zero variance, 5^, does not alter 
the prediction of the best-fit value of /3o, although it nat- 
urally does lead to a broadening of the areas of parameter 
space included at the 2u level. There is relatively little 
broadening of the region of parameter values included at 
the 1(7 level in the /3o = /?* plane, however, which is no 
doubt a consequence of the steepness of the go contours 
in Fig. 0(a) in the area corresponding to region II. 

In Figs, g and ^ we display the joint credible regions 
for M and /, for two slices through the 3-dimensional 
{M,f,l3o) parameter space for Wi = 1.5: (a) the best-fit 
case /3o = /?*; and (b) /3o = 0. Analogously, the best-fit 
case is shown in Fig. [ill for Wi = 0.2. We see from Fig. 

that once the likelihood is normalized relative to (3* 
no regions remain in the (3 — parameter plane which 
are admitted at the Icr level when b ~ 0. Furthermore, 
even when a non-zero standard deviation, 6, is included, 
region II of the (M, /) parameter plane is favoured at the 
Icr level, in contrast to Fig. |[ 

We have also undertaken an analysis of the models 
with (3o — but variable 6, similarly to the study of Ref. 
p5[ . In that case, we once again find that region II of 
the (il/, /) parameter plane is admitted at the la level if 
b = 0.25 or b = 0.5. The dependence of the value of 



Q = -21n(/:Aiy) = - 21n 



i^/K sa 



on the value of b is displayed in Fig. |l^, for /3o = as 
compared with the best-fit case (3q = P* . The quantity Q 
is analagous to the chi-square statistic of the maximum 
likelihood method. We see that the /3o = models favour 
a non-zero value of 5 ~ 0.36 by a very small margin as 
compared to the b — case. For b — 0.36 the points of 
greatest likelihood lie mainly in region II, in contrast to 
the 5 = case in Fig. @. 
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FIG. 12. Variation of the least value of Q = — 21n(£A!/) 
as a function of the standard deviation b of the prior distri- 
bution for for values /3o = and jSo = {3* . 

Varying the inital condition Wi does not appear to af- 
fect the best-fit value of /3o significantly. For Wi — 0.2 
(cf.,Fig. ^l]), for example, the best-fit value was (3* — 
0.435, a difference of 5% from the Wi — 1.5 case. Fur- 
thermore, the numerical value of the least value of Q was 
only 0.2% greater in the Wi ~ 0.2 case. We have not 
attempted to find a best- fit value for Wi. 



V. CONSTRAINTS FROM LENSING STATISTICS 

Gravitational lensing of distant light sources due to the 
accumulation of matter along the line of sight provide an- 
other relatively sensitive constraint on the cosmological 
models of interest. For cosmology the situation of most 
interest is the lensing of high luminosity quasars by in- 
tervening galaxies. The abundance of multiply imaged 
quasars and the observed separation of the images to the 
source puts constraint on the luminosity-redshift relation 
and hence the model parameters. Basically, if the volume 
of space to a given redshift is larger then on average one 
can expect more lensing events. This leads to a statisti- 
cal test, which has been used to put bounds on A (2^-|2^] 
and to test properties of some decaying A or quintessence 
models 

Gravitational lensing statistics are useful since they 
provide a tests which potentially provides opposing con- 
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(b) = 0.2 M(lO-V/^eV) 

FIG. 13. Confidence limits from gravitational lensing statistics: (a) Wi = 1.5; (b) Wi = 0.2. Parameter values excluded at 
the 95.4% level are darkly shaded, while those excluded at the 68.3% level axe lightly shaded. For reference, contours of fi^o 
and Hoto are superposed as dashed and dotted lines respectively. 
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straints to those obtained from supernovae magnitude- 
redshift tests. In particular, in the case of models with 
a vacuum energy provided by a cosmological constant, 
the high redshift supernovae have been interpreted as 
favouring relatively large values of J^a - Perlmutter et al. 
^ give a value of JIa ~ 0.72 at la level - whereas the 
gravitational lensing data leads to upper bounds on ^Ia: 
Kochanek |^ quotes fl\ < 0.66 at the 2a level. A com- 
bined likelihood analysis has been performed by various 
authors ||,|^J30|. 

Gravitational lensing constraints on the PNGB mod- 
els have been very recently given by Frieman and Waga 
for Wi = 1.5. However, Frieman and Waga consid- 
ered a more retricted range of parameter space, M < 
O.OObh eV , since they did not consider the possibility 
of source evolution in the case of the type la supernovae 
and therefore took values of M > 0.005/i eV to be ruled 
out. We wish to extend the range of M for the grav- 
itational lensing statistcs to consider parameter values 
corresponding to region II in the supernovae constraint 
graphs, Figs. |o|, so as to compare the constraints from 
different tests. 

We have thus simply followed the calculation described 
by Waga and Miceli who performed a statistical 

lensing analysis of optical sources described earlier by 
Kochanek [g5[. They used a total of 862 (z > 1) high lu- 
minosity quasars plus 5 lenses from seven major optical 
surveys ||38| . (Another alternative not considered here is 
to analyse data from radio surveys - see, e.g., p6| , p7t .) 
Undertaking a similar analysis for the increased param- 
eter range, we arrive at Fig. |l^, which shows the 68.3% 
and 95.4% joint credible regions for M and /, for two val- 
ues of Wi. We refer the reader to Refs. p5| , ^ for details 
of the calculation. The only regions of parameter space 
excluded at the 2a level turn out to be areas of parame- 
ter space for which the deceleration is presently negative 
(cf.. Fig. 0), with the scalar field still commencing its first 
oscillation at the present epoch. 

VI. DISCUSSION 

Let us now consider the overall implications of the con- 
straints observed above. 

Firstly, since empirical models with source evolution do 
appear to fit the data somewhat better, it would appear 
that we do have weak evidence for an underlying evo- 
lution of the peak luminosity of the type la supernovae 
sources, at least in the context of the PNGB quintessence 
models. It might be interesting to compare the case of 
other quintessence models, or the case of a cosmological 
constant. However, the PNGB model is qualitatively dif- 
ferent from such models since its final state corresponds 
to one in which the ultimate destiny of the universe is to 
expand at the same rate as a spatially flat Friedmann- 
Robertson-Walker model, rather than to undergo an ac- 
celerated expansion. This is of course precisely why we 



chose the PNGB models as the basis of our investigation, 
rather than models in which a late-time accelerated ex- 
pansion had been built in by hand. If we wish to test 
the hypothesis that the faintness of the type la super- 
novae is at least partly due to an intrinsic variation of 
their peak luminosities - which is a very real possibility 
in view of the results of ||^ - then a quintessence model 
which possesses a variety of possibilities for the present- 
day variation of the scale factor is probably the best type 
of model to investigate. 

If only supernovae luminosity distances, (cf.. Fig. 
and gravitational lensing statistics, (cf.. Fig. are 
compared then we see that there is a remarkable con- 
cordance between the two tests - region II of Fig. I com- 
cides with a region included at even the la level in Fig. 

P. This is perhaps not surprising, since in view of Fig. 
region II corresponds to parameter values for which 
the present day universe has already undergone almost 
one complete oscillation of the scalar field about the fi- 
nal critical point C2 of Fig. |^. It is thus already well on 
the way towards its asymptotic behaviour, which closely 
resembles that of a standard spatially flat Friedmann- 
Robertson- Walker model. 

Due to the oscillatory behaviour, parameter values in 
region II correspond to models in which there has been 
a recent cosmological acceleration, (e.g., at Z'--^0.2), but 
with a q{z) which changes sign three times over the larger 
range of redshifts, < 2: < 4, in the quasar lensing sam- 
ple, and therefore differing significantly from Friedmann- 
Lemaitre models over this larger redshift range. Extend- 
ing the SNe la sample to include objects at redshifts 
0.1 < z < 0.4 and z > 0.85 in substantial numbers would 
greatly improve the ability to decide between models in 
regions I and II. 

There is some cause for concern, however, if we con- 
sider the favoured values of O^o and -ffo^o- In the case 
of the models with empirical evolution of the supernovae 
sources, we always found that the best-fit parameter val- 
ues occured at the il^o 1 boundary of the {M, f) pa- 
rameter space. The overwhelming evidence of many as- 
tronomical observations over the past two decades |3|] 
would tend to indicate that 0^0 — 0.2 ± 0.1, indicating 
that a vacuum energy fraction of fi^o^ 0.7-0.8 is desir- 
able, and fl^o < 0.9 in any case. Although parameter 
values with Q^q < 0.8 certainly fall within both the 2a 
and la portions of region II of Fig. ^ for all values of 
b, there are potentially serious problems if we wish to 
simultaneously obtain large values of HqIq. In view of 
recent estimates of the ages of globular clusters a 
lower bound of 12 Gyr for the age of the Universe ap- 
pears to currently indicated. With h ~ 0.65 this would 
require HqIq > 0.8. For Wi = 1.5, parameter values with 
Hoto > 0.8 coincide with values fi^o ^ 0.9 in region II, 
which is phenomenologically problematic. 

The tension between the values of Q^q and Hoto is 
somewhat mitigated for lower values of Wi. For Wi = 0.2, 
for example, we see from Figs. ^ and |l|(b) that the 
^4>o = 0.7 and Hoto meet in Region II, and there is a 
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small region of parameters there with 0.7 < fi^o < 0.9 
and Hoto > 0.8, which is also consistent with the other 
cosmological tests. 

Even if the supernovae sources undergo evolution it is 
clear that parameter values in region I of I which are 
favoured in the absence of evolution of peak Snia lumi- 
nosities, are still included at the 2(7 level in the models 
with evolution, in view of Fig. ^. 

Perhaps the most significant aspect of our results is 
the fact that the mere introduction of an additional dis- 
persion, b > 0.17, in the peak luminosities while leaving 
their mean value fixed (cf. , Fig. |l^) , gives rise to a change 
in the best-fit region of parameter space from region I to 
region II. (See pdj for further details.) One would imag- 
ine that an increased dispersion is likely to be a feature 
of many models of source evolution, even if evolution- 
ary effects are of secondary importance. Thus even if 
the empirical models with non-zero /3 are somewhat ar- 
tificial, more sophisticated scenarios could well lead to 
similar changes in regard to the fitting of cosmological 
parameters in the PNGB model. 

Much tighter bounds on the parameter space of 
quintessence models, including the present model, will 
be obtained over the next decade as more supernovae 
data is collected. What we wish to emphasize, however, 
is that an effective vacuum energy which is cosmologi- 
cally significant at the present epoch should not simply 
be thought of in terms of a "cosmic acceleration" . A dy- 
namical vacuum energy with a varying effective equation 
of state allows for many possibilities for the evolution 
of the universe, and overly restrictive assumptions, such 
as equating quintessence to models with a late period of 
continuous cosmological acceleration, should be avoided. 
If detailed astrophysical modeling of type la supernovae 
explosions ultimately shows that the dimness of distant 
supernova events is largely due to evolutionary effects, 
it does not spell the end for cosmologies with dynamical 
scalar fields. 
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